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Abstract. We study the effects of the allocation of distributed gener- 
ation on the resilience of power grids. We find that an unconstrained 
allocation and growth of the distributed generation can drive a power 
i-£h ' grid beyond its design parameters. In order to overcome such a problem, 

we propose a topological algorithm derived from the field of Complex 
Networks to allocate distributed generation sources in an existing power 
grid. 
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Distributed Generation from renewable sources is having a deep impact on our 
power grids. The difficult task of integrating the stochastic and often volatile 
renewable sources into a the grid designed with a power-on-demand paradigm 
could perhaps solved leveraging on distributed storage [5]; nevertheless, massive 
and economic power storage is not yet readily available. In the meanwhile, power 
grids are nowadays required to be robust and smart, i.e. systems able to maintain, 
under normal or perturbed conditions, the frequency and amplitude variations 
of the supplied voltage into a defined range and to provide fast restoration after 
faults. Therefore, many studies have concentrated on the dynamic behaviour of 
Smart Grids to understand how to ensure stability and avoid loss of synchro- 
nization during typical events like the interconnection of distributed generation. 
The large number of elements present into real grids calls for simplifications 
like the mapping among the classic swing equations [53] and Kuramoto models 
[17118116] that allows to study numerically or analytically the synchronization 
and the transient stability of large power networks. 

Even simple models [TH] akin to the DC power flow model [25] show that 
the network topology can dynamically induce a complex size probability distri- 
butions of blackouts (power-law distributed) , both when the system is operated 
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near its limits [TT] or when the system is subject to erratic disturbances [22] . 
New realistic metrics to assess the robustness of the electric power grid with 
respect to the cascading failures [26] are therefore needed. 

Smart grids are going to insist on pre-existing networks designed for differ- 
ent purposes and tailored on different paradigms and new kind of failures are 
possible: therefore a careful transition is needed. One possible approach could 
be the use of advanced metering infrastructure (AMI) not only for implement- 
ing providers and customers services, but also to detect and forecast failures; 
nevertheless an ill-designed network will never be efficient. 

Our approach will not concentrate on the instabilities but will focus instead 
on the condition under which, in presence of distributed generation, the system 
can either be operated or controlled back within its design parameters, i.e. it is 
resilient. It is akin in spirit to the approach of [13] . that by applying DC power 
flow analysis to a system with a stochastic distribution of demands, aims to 
understand and prevent failures by identifying the most relevant load configu- 
rations on the feasibility boundary between the normal and problematic regions 
of grid operation. 

To model power grids, we will use the more computational intensive AC 
power flow algorithms as, although DC flows are on average wrong by a few 
percent [24], error outliers could distort our analysis. 

To model distributed renewable sources, we will introduced a skewed prob- 
ability distribution of load demands representing a crude model of reality that 
ignores the effects like the correlations (due for examples to weather conditions) 
between different consumers or distributed producers . 



2 Methods 

2.1 AC Power Flow 

The AC power flow is described by a system of non-linear equations that allow 
to obtain complete voltage angle and magnitude information for each bus in a 
power system for specified loads [TH] ■ A bus of the system is either classified as 
Load Bus if there are no generators connected or as a Generator Bus if one or 
more generators are connected. It is assumed that the real power PD and the 
reactive power QD at each Load Bus are given, while for Generator Buses the 
real generated power PG and the voltage magnitude \V\ are given. A particular 
Generator Bus, called the Slack Bus, is assumed as a reference and its voltage 
magnitude \V\ and voltage phase are fixed. The branches of the electrical 
system are described by the bus admittance matrix Y with complex elements 
YijS. 

The power balance equations can be written for real and reactive power for 
each bus. The real power balance equation is: 

N 

= -P t +Yl W \V k \ (G tk cose ik + B lk sin**) 

k=l 
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where N is the number of buses, P, is the net real power injected at the i 
bus , Gik is the real part and Bn~ is the imaginary part of the element Yjj and 
Oik is the difference in voltage angle between the i th and k th buses. The reactive 
power balance equation is: 

N 

= + \Vi\ \V k \ (G tk sm6. lk - B lk cos9 tk ) 
k=i 

where Qi is the net reactive power injected at the i th bus. 

Real and reactive power flow on each branch as well as generator reactive 
power; the output can be analytically determined but due to the non-linear 
character of the system numerical methods are employed to obtain a solution. 
To solve such equations, we employ Pylon [2J, a port of MATPOWER jj] to the 
Python programming language. 

A requirement for the stability of the load and generation requirements is 
the condition that all branches and buses operate within their physical feasibil- 
ity parameters; going beyond such parameters can trigger cascades of failures 
eventually leading to black outs [21] , 

In the present paper a topological investigation on the power grid has been 
developed in order to evaluate the effects of distributed generation on the voltage 
and power quality. Hence, a steady state analysis has been carried out and the 
transient phenomena connected to the power flow control have been neglected. 
Under this hypothesis the frequency variation connected to power flow control 
has been considered stabilized and the system has been considered characterized 
by a constant steady state supply voltage frequency. Therefore, if all the nodes 
are near their nominal voltage, it is much easier to control the system and to 
avoid reaching infeasible levels of power flow. Consequently, to measure the ef- 
fects of power quality of a power grid under distributed generation we measure 
the fraction F of load buses whose tension goes beyond ±5% of its nominal 
voltage. Notice that real networks are often operated with some of the buses 
beyond such parameters so that (especially for large networks) it is expected to 
be F ^ under operating conditions. The maximum of the resilience for a power 
grid (intended as the capability of restoring full feasible flows) is expected to be 
for F = 0. 

2.2 Distributed Generation and Skew-normal distribution 

We will consider distributed generation due to erratic renewable sources like 
sun and wind; therefore, we will model the effects of "green generators" on a 
power grid as a stochastic variation the power requested by load buses. Load 
buses with a green generator will henceforth called green buses. We will consider 
the location of green buses to be random; the fraction p of green buses will 
characterize the penetration of the distributed generation in a grid. 

If the power dispatched by distributed generation is high enough, loads can 
eventually become negative: this effect can be related to the efficiency of green 
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generators. We model such an effect by considering the load on green buses 
described by the skew- normal distribution [3] , a pseudo- normal distribution with 
a non-zero skewness: 

/ (x, a) = 2cj) (x) <P (ax) 
where a is a real parameter and 

<t> (x) = exp (~x 2 /2) /v^tt <2> (ax) = f <j> (t) dt 

The parameter a will characterize the level of the distributed generation: to 
positive a correspond loads positive on average, while for negative a green nodes 
will tend to dispatch power. 

Our model grids will therefore consist of three kind of buses: Nq generators 
(fixed voltage), Ni pure loads (fixed power consumption) and N g green buses 
(stochastic power consumption) with Nq + Ni + N g = N the total number of 
buses and N g + Ni — Nl the number of load nodes. The fraction p = N 9 /Nl 
measures the penetration of renewable sources in the grid. 



2.3 Complex Networks and Page Rank 

The topology of a power grid can be represented as a directed graph G = (V,E), 
where to the i-th bus corresponds the nodes n,; of the set V and to the fc-th 
branch from the i-th to the j-th bus corresponds the edge = (i,j) of the 
set E. In Power System engineering, it is custom to associate to the graph G 
representing a power networks its incidence matrix B whose elements are 

f life k = (i,.)eE 
B lk = I -lif e k = (_, i) e E 
[ otherwise 

. An alternative representation of the graph much more used in other scientific 
fields is its adjacency matrix A whose element are 

13 \0 otherwise 

While Graph Theory has an old tradition since Euler's venerable problem 
on Koenigsberg bridges [5] , Complex Networks is the new field investigating the 
emergent properties of large graphs. An important characteristic of the nodes 
of a complex network is their centrality, i.e. their relative importance respect 
to the other nodes of the graph [TU]. An important centrality measure is Page 
Rank, the algorithm introduced Brin and Page [7] to rank web pages that is at 
the hearth of the Google search engine. The Page Rank rj of the i-th node is the 
solution of the linear system 

N d° 
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where N is the number of buses (nodes), d° = Aij is the number of outgoing 

i 

links (out-degree) and p = 0.85 is the Page Rank damping factor. In study- 
ing power grids, we will employ Page Rank as it is strictly related to several 
invariants occurring in the study of random walks and electrical networks [14] . 

3 Results 

3.1 Effects of distributed generation 

We have investigated the effects of our null model of distributed generation on 
the 2383 bus power grid of Poland, 1999. Starting from the unperturbed net- 
work, we have found an initial fraction Fq = 1.6% of load buses beyond their 
nominal tension. We have therefore varied the penetration p at fixed distributed 
generation level a's; results are shown in Fig. Q] 
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Fig. 1. Effects of the penetration p of distributed generation on the resilience of the 
Polish power grid at different values a of the green generators. Notice that for a = 
renewable sources satisfy on average the load requested by the network, while for a < 
there is a surplus of renewable energy. Lower values of the fraction F of buses operating 
near their nominal tension correspond to a higher resiliency. Notice that the penetration 
of distributed generation initially enhances resiliency. At higher values of p, resilience 
worsens; in particular, it is severely impaired if distributed generation produces on 
average more energy than the normal load requests (a = — 1). It is therefore advisable 
to keep levels of renewable energy production below the normal load request (a = 1). 

We find that the behaviour of the fraction F of buses operating near their 
nominal tension does not follow a monotonic behaviour. Initially (low values of 




6 Antonio Scala et al. 



p), the penetration of distributed generation enhances resiliency (i.e. decreases 
F). At higher values of p, F grows and resilience worsens. Such an effect is 
particularly severe if green nodes introduce a surplus (a < 0) of power respect 
to the normal (p = 0) operating load requests. On the other hand, keeping 
the levels of renewable energy production below (a > 0) the normal load request 
delays the point beyond which the penetration of distributed generation worsens 
the resiliency. 

Notice that when distributed generation is ancillary (a > 0) and not predom- 
inant in the power supplied of the network, full penetration {p = 1) of renewable 
sources lead to more stable state than the initial (p = 0) one. 

3.2 Targeted distributed generation 

Beside their natural application to web crawling, the Page Rank algorithm can 
be applied to find local partitions of a network that optimize conductance [3] . We 
therefore investigate what happens in a power network if distributed generation 
is introduced with a policy that accounts for the pagerank of load nodes. In other 
words, for a level of penetration p, we choose the first n g = pN^ load nodes in 
decreasing pagerank order to become green nodes. The effects of such a choice 
are shown compared to the random penetration policy in Fig. [5] 

We find that, for low penetration levels, the pagerank policy reduces the 
number of nodes operating beyond their nominal tension both for positive and 
for negative a's. Again, the excess of power production (a < 0) comparatively 
reduces the resilience of the network. 

Preliminary results show that Page Rank is the best behaved among cen- 
tralities in enhancing power grid resilience; such study will be the subject of a 
future publication. 

4 Discussion 

We have introduced a model base on the AC power flow equation that allows to 
account for the presence of erratic renewable sources distributed on a power grid 
and for their efficiency. By defining the resilience of the grid as a quantity related 
to the possibility of controlling the power flow via voltage adjustments (hence 
returning within the operating bounds of its components), we have studied the 
penetration of distributed generation on a realistic power grid. 

We have found that while the introduction of few "green" generators in gen- 
eral enhances the resilience of the network by decreasing the number of nodes 
operating beyond their nominal voltage, a further increase of renewable sources 
could decrease the power quality of the grid. Anyhow, if distributed generation 
is ancillary and not predominant in the power supplied of the network, the grid 
at full penetration (p — 1) of renewable sources is in a more stable state than 
the starting grid (p = 0). 

Our finding that a surplus of production from renewable sources is also a 
source of additional instabilities is perhaps to be expected in general for networks 
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Fig. 2. Comparison between random placement (filled symbols) and page- rank place- 
ment (empty symbols) of green generators in the Polish grid, both for surplus produc- 
tion of renewable energy (upper panel, a = — 1) and for levels of renewable energy 
production below the normal load request (lower panel, a — 1). The page- rank place- 
ment of renewable sources allows to attain lower values of the fraction F of buses 
operating near their nominal tension (and hence a higher resiliency) at lower values of 
the penetration p. The best case is realized for levels of renewable energy production 
below the normal load request, where a plateau to low values of F is quickly attained. 
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that have been designed to dispatch power from their generators to their loads 
and not to produce energy "locally". While we have found that in an isolated 
grid instability possible increases with the penetration, what happens when more 
grids are linked together is an open subject. Power grids are typical complex 
infrastructural systems; therefore they can exhibit emergent characteristics when 
they interact with each other, modifying the risk of failure in the individual 
systems |12) . As an example, the increase in infrastructural interdependencies 
could either mitigate [8 or increase [20 9 the risk of a system failure. 

Finally, we find that a policy of choosing the sites where to introduce re- 
newable sources according to Page Rank allows to increase the resilience with 
a minimal amount of green buses. Such policy does not take into account other 
factors and should therefore be integrated in a multi-objective optimization to 
consider the environmental, economical and social constraints. 
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